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It is well known that the discrete Sierpinski triangle can be defined as the nonzero 
residues modulo 2 of Pascal's triangle, and that from this definition one can easily 
construct a tileset with which the discrete Sierpinski triangle self-assembles in Winfree's 
tile assembly model. In this paper we introduce an infinite class of discrete self-similar 
fractals that are defined by the residues modulo a prime p of the entries in a two- 
dimensional matrix obtained from a simple recursive equation. We prove that every 
fractal in this class self-assembles using a uniformly constructed tileset. As a special 
case we show that the discrete Sierpinski carpet self- assembles using a set of 30 tiles. 



1 Introduction 

A model for self-assembly is a computing paradigm in which many small components interact 
locally, without external direction, to assemble themselves into a larger structure. Wang 
[HI E] ■ first investigated the self-assembly of patterns in the plane from a finite set of square 
tiles. In Wang's model, a tile is a square with a label on each edge that determines which 
other tiles in the set can lie adjacent to it in the final structure. The Tile Assembly Model 
of Winfree pi)] , later revised by Rothemund and Winfree [HI [5] , refines the Wang model to 
provide an abstraction for the physical self-assembly of DNA molecules. 

We introduce some formal notation for the Tile Assembly Model in the next section. Briefly, 
a tile is a square with a label on each edge, which we represent as a string, but in addition 
each edge has a integer bonding strength of 0, 1, or 2, represented in Figure [T] by a dashed 
line, a solid line, or a double line, respectively. A tile may also have a label in the center 
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for informational purposes. Tiles are assumed not to rotate. Two tiles can potentially lie 
adjacent to each other only if the adjacent edges have the same label and the same bonding 
strength. Intuitively, the bonding strength and edge label model the bonding strength and 
"sticky ends" of a specially constructed DNA molecule, as illustrated in Figure [T] 

A tile system is assumed to start with an infinite supply of a finite number of tile types. A 
set of initial tiles, the seed assembly, is placed in the discrete plane. Self-assembly proceeds 
nondeterministically as new tiles bond to the existing assembly. The ability of tiles to bond 
is controlled by a system parameter called the temperature. In this paper we are concerned 
with temperature 2 systems, which means that a tile may bond to an existing assembly only 
if the sum of the bonding strengths of the edges in the tile that abut the assembly is at least 
2. 

In this process, tiles may cooperate to create a planar structure, or (by appropriate inter- 
pretation of the labels) perform a computation. Winfree [10] and others [6j El [2j [TJ have 
shown that such systems can perform computations such as counting and addition, and that 
in fact the model is universal: given an arbitrary Turing machine, there is a tile set for 
which each row of the resulting assembly is the result of a computation step of the Turing 
machine. It is also possible to use a finite tile set to generate infinite planar structures such 
as discrete fractals. The latter was made famous when Papadakis, Rothemund and Winfree 
[I] performed an experiment in which actual DNA molecules were used to self-assemble a 
portion of the discrete Sierpinski triangle. 

The discrete Sierpinski triangle has been used extensively as a test structure for in DNA 
self-assembly [TU]. One reason for this is that it self-assembles using a simple set of only 7 
tiles [ID] • More generally, however, fractal structures are of interest because "Structures that 
self-assemble in naturally occurring biological systems are often fractals of low dimension, by 
which we mean that they are usefully modeled as fractals and that their fractal dimensions 
are less than the dimension of the space or surface that they occupy. The advantages of 
such fractal geometries for materials transport, heat exchange, information processing, and 
robustness imply that structures engineered by nanoscale self-assembly in the near future 



2 



will also often be fractals of low dimension." [2] It is then natural to ask what other discrete 
fractals, other than the ubiquitous Sierpinski triangle, can self-assemble with a relatively 
small set of tile types in this model. 

In this paper we introduce an infinite class of self-similar discrete fractals, all of which self- 
assemble in Winfree's model. The class includes, as special cases, the standard Sierpinski 
triangle and Sierpinski carpet. All the fractals in this class exhibit a strong self-similarity 
property that we call numerical self- similarity. Each fractal is defined in terms of an infinite 
integer matrix M whose entries are residues modulo a given prime q. A fractal S (as a subset 
of the first quadrant of integer plane) can then be defined as the set of points for which 
M[i,j] is not congruent to zero, modulo q. The usual notion of self-similarity for a set S 
means that there is an integer p such that for any s,t < p, the set of points consisting of the 
p k by p k square whose lower left corner is at sp k ,tp k is either empty, if (s,t) S, or is an 
exact copy of the p k by p k square wiht lower left corner at the origin, if (s, t) is in S. 

Numerical self-similarity means further that the entries of the p k by p k submatrix of M with 
lower left corner at M[sp k ,tp k ] are always related to those of the p k by p k submatrix at the 
origin by a factor of M[s, t], that is, 

M[sp k + i,tp k + j] = M[i,j]M[s,t] 

for i,j < p k . One consequence of our result is that there exists a simple recursively defined 
matrix M that defines the discrete Sierpinski carpet, using p = 3. Figure [T] is an illustration 
of the discrete Sierpinski carpet as a subset of the plane, and Figure [I] is a depiction of the 
mod 3 residues of the matrix M that defines it. It follows from the simple definition of the 



matrix M and Theorem 3J3 that the Sierpinski carpet self-assembles in the Tile Assembly 
Model. 

The next section introduces some definitions and notation for the Tile Assembly Model 
described above. Section [3] reviews some known results on the relationship between the 
Sierpinski triangle and Pascal's triangle and gives a uniform construction of tilesets for 
recursively defined matrices. In Section [4] we define the matrices from which we obtain 
discrete fractals and prove the main result on numerical self-similarity of such matrices. 
Section [5] contains some concluding remarks and open problems. 



2 Preliminaries 



In this section we introduce some notation and terminology associated with the Tile Assembly 
Model described in the introduction. The description given here should be sufficient for our 
purposes; for formal details see 0, El EE]- 

We work in the discrete Euclidean plane. Let U2 denote the set of unit vectors, denoted 
as cardinal directions N = (0,1), S = (0,-1), E = (1,0), and W = (—1,0). Let S be a 
finite alphabet. A tile type t is a a pair of functions (col t ,str t ) where cob : U2 — > S* and 
col t : U2 — > N. That is, a tile type associates a color col t (u) G £* and a strength str t (u) with 
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Figure 2: First three stages of the Sierpinski carpet. 
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Figure 3: First three stages of the numerically self-similar Sierpinski carpet. 
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each of the four sides of a unit square, where the side is indicated by the unit vector u. We 
also assume that there is a label associated with each tile type by a function m : T — > C, 
where £ is a finite alphabet (which in the examples of interest will be a set {0, 1, ... ,p — 1} 
for some prime p) . 

Let T denote a finite set of tile types and let r e N be a fixed parameter, called the 
temperature (which in the present paper is always 2). In general we assume that there is 
an infinite supply of tiles for each type t e T. A tile may be positioned, but not rotated, 
in the discrete plane. As described in the introduction, two adjacent tiles may bond if the 
abutting edges have matching color and matching strength s; the strength of the bond is s. 
More generally, an assembly is a partial assignment a : Z — > T of tile types to locations in 
the plane in which each tile is bonded to its neighbors with a total strength of at least r, 
and such that the assembly cannot be separated into smaller assemblies without breaking a 
set of bonds having a total strength of at least r. 

The process of self-assembly begins with a given seed assembly o and proceeds nondetermin- 
istically by extending the domain of the assembly, where a new tile may extend an assembly 
at position (x,y) if all edges abutting those of existing adjacent tiles have matching colors 
and matching strengths and if the sum of the strengths for the abutting edges is at least r. 
An assembly is terminal if it cannot be extended. 

A tile assembly system (TAS) is a triple (T, a, r) where T is a finite set of tile types, a is 
the seed assembly, and r is the temperature. A TAS is definitive if it has a unique terminal 
assembly. 

We define self-assembly for first for matrices with values in an arbitrary alphabet, and then 
for subsets of the discrete plane. 

Definition 2.1 Fix an alphabet C. For any matrix M with values in C and any subset £ 
of C, the set S C Z 2 determined by (M,C) is the set of points (x,y) such that M[x,y] is 
defined and M[x,y] G C . 

Definition 2.2 (a) Let C be a finite alphabet and M a matrix, possibly infinite, with values 
in C M self-assembles if there exists a definitive TAS (T, a, r) with terminal assembly 
a, and a labeling m : T — > C, such that a(x,y) is defined if and only if M[x,y] is 
defined and for all (x,y) in the domain of M, m(a(x,y)) = M[x,y]. 

(b) Let S C Z 2 . S self-assembles if there exists a finite alphabet C, a matrix M with values 
in C, and a subset C of C such that M self- assembles and S is the subset determined 
by (M,£>). 
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3 Tiling the Sierpinski triangle 



In this section we formally introduce the notion of numerical self-similarity mentioned in 
the introduction and review some known results regarding self-assembly of the Sierpinski 
triangle. We then present a uniform construction of tile assembly systems for recursively 
defined matrices. 

The following definition generalizes the usual definition of a discrete self-similar fractal. 



Definition 3.1 Let p, q > 2. Let M : N 2 — > N. M is numerically p-self-similar modulo q if 

for all < s,t < p, for all k > and for all i,j < p k , 

M[sp k + i,tp k +j] = M[s,t\- M[i,j] (1) 

modulo q. If M is a finite matrix, M is numerically p-self- similar modulo q if holds 
wherever M is defined. 



Suppose a matrix M defined for all i, j > has values in {0, 1, . . . , q — 1}. It is not difficult 
to see that if M is numerically p-self-similar modulo q and S C N 2 is the set defined by 



(M, {1,2, ... q — 1}) in the sense of Definition 2.1 (i.e., (x, y) E S if and only if M[x, y] ^ 



modulo q), then S is a discrete self-similar fractal in the usual sense. 

We define Pascal's triangle to be the matrix P of integers defined for all i,j > by 

P[i,0] = 1 
P[0J] = 1 

P[i,j] = P[i,j-1]+P\j,i-1] forz,j>0 



Thus P[i,j] — |^ . j, where the reverse diagonal i + j = k corresponds to the usual kth 

row of Pascal's triangle rendered horizontally. The relationhip between Pascal's triangle and 
the Sierpinski triangle is well known; if M is a matrix containing the residues modulo 2 of 
P, then the set determined by M and {1} is the Sierpinski triangle. This fact, it turns out, 
is a special case of a more general result. 



Theorem 3.2 P is numerically p-self- similar modulo p for any prime p. 



Theorem 3.2, in turn, is a special case of our Theorem 4.2 It has been known in various 
forms for some time; for example, there is a proof in [3], and a form of it is implicit in 
Winfree's 7-tile system in which the Sierpinski triangle self-assembles |10j . 

In the remainder of this section we present a uniform construction of a tile assembly system 
for any recursively defined matrix, in the following sense. Given an matrix M, if there is a 
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positive integer n such that M[x,y] is determined by a finite function of the entries M[x',y'] 
with x — n < x' < x and y — n < y' < y (excluding M[x,y] itself, of course), then there 
is a tile assembly system in which M self-assembles. This claim is formalized in the next 
theorem. 

Let M : Z 2 — > C be an infinite two-dimensional matrix whose entries come from some finite 
alphabet C For x, y G Z, define 

R M (x,y) = [fi,f 2 , ■ ■ . ,f n _i] T , 

where each row vector r*j is defined as 

n = (M[x-i,y-n + l],...,M[x-i,y-l],M[x-i,y]) 

for 1 < i < n, and define 

f M (x,y) = (M[x,y-n + l],...,M[x,y-l}). 

That is, R,M(x,y) is the n-lxn submatrix whose upper right corner is at (x — l,y), and 
rM{x,y) is the vector consisting of the n — 1 elements of row x directly to the left of (x,y). 
Thus the pair (rM{x,y), Rm{x,u)) contains the entries of the n x n submatrix of M whose 
upper-right corner is at (x, y), excluding M[x, y] itself. Taking the entries in row-major order 
we can identify (f M (x,y), RM(x,y)) with an element of C n -1 . 

Theorem 3.3 Let C be any finite alphabet and let _L be a symbol not in C Let C± = £U{_L}. 

Let f be any function 

f : C ± n ' 2 - 1 -> C. 

Define the matrix M : Z 2 — > C± by 

M[x,y] = ±ifx<0ory<0 

M[x,y] = f(f M (x,y),R M {x,y)) otherwise. 

Then there is a definitive tile assembly system T = (T, a, 2), and a labeling m : T — > £ 7 sttc/j 
i/iai m i/ie unique terminal assembly a of T , 

m(a(x,y)) = M[x,y] 

for all x, y > 0. 

The construction of the tileset is as follows. For each element (r, R), we define a tile type t. 
Let 

f = (r ,r 1? . . . ,r„_ 2 ), 
# = [ri,f 2 , . . . ,f„_i] T , 
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and define 

b = f(f,R), 

r> = (ri,r 2 , . . . ,r n - 2 ,b), 

r = (r ,ri, . . . ,r n - 2 ,b), and 

R' = [r ,r 1 ,r 2 ,...,r n _ 2 ] T . 

Then the color of t — t(f, R) is defined by 



co\ t (W) 


= f 


co\ t (S) 


= R 


co\ t (E) 




co\ t (N) 


= R' 



and tile t is labeled as m(t) = b. 

Finally, the strength of t is str t (u) = 1 for all unit vectors u except for the following cases: 

Seed tile: If = (±, . . . , _L) for 1 < % < n and r = (±, . . . , _L), then str t (JV) = 2 and 
str t (£7) = 2. 

Row 0: If fi = (±, . . . , _L) for 1 < % < n but f ^ (_L,...,_L), then str t (W") = 2 and 
str t (£) = 2. 

Column 0: If r*j = (_L, _L, . . . , _L, r), where r ^ _L, for 1 < i < n and r = (±, . . . , _L), then 
str t (5) = 2 and str t (7V) = 2. 

A typical tile is shown below. 

R! 

r b r' 
R 

Define the seed assembly for T by er(0, 0) = s, where s is the unique seed tile described 
above. 

The proof of correctness for the above construction is straightforward and is omitted. 
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4 Tiling the Sierpinski carpet 



In this section we introduce a class of matrices defined by a simple recursion and prove 
our main technical result, which states simply that matrices in this class are numerically 
self-similar. As a consequence, we are able to conclude in Corollary |4.3| that the Sierpinski 
carpet self-assembles. 



Definition 4.1 Let a, b, c> 0. Let M be defined for all > by 

M[0,0] = 1 

M[Q,j] = a j forj>0 

M[i, 0] = e for i > 

M[i,j] = aM[i, j - 1] + bM\i - 1, j - 1] + cM[i - for i,j > 



(2) 



When a = b = c = 1, the entries of M are known as the Delannoy numbers (see the 
Encyclopedia of Integer Sequences A001850 [7] ). 

Theorem 4.2 M is numerically p- self- similar modulo p. 



Before moving to the proof of Theorem 4.2 we briefly discuss the following consequence. 



Corollary 4.3 The Sierpinski carpet self-assembles. 



Proof. Let a = b = c = 1 and p = 3 in Definition 4A Note we can assume that the entries 
of M are the residues modulo 3. Let 5cNbe the set of points (x, y) such that M[x, y] ^ 0. 



Since S is self-similar, evidently S is the Sierpinski carpet. By Theorem 3.3, S self-assembles. 
□ 



Note that the tileset constructed by Theorem 3.3 is not optimal; in particular most of the 
tiles involving the symbol _L are unused. It is not difficult to show that a set of 30 tiles is 
sufficient. 

Define a tileset T to consist of the three tiles 



:m) 
i i 



i i i 



i i 
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plus 27 tiles of the form 




for x,y,z G {0, 1, 2}, where w = x + y + z (mod 3). 



We next proceed to prove some preliminary results supporting Theorem ^2 The following 
is a very slight generalization of a well-known combinatorial interpretation of the Delannoy 
numbers. We give a brief argument here for completeness. Note that a similar expression 
holds when % > j. 



Definition 4.4 For i < j let 



k=0 



j \ ( j+i-k 
i — k 



(3) 



Claim 4.5 For i < j , 



M[i,j] = D{i,j). 



(4) 



Proof. Define a move as a pair of coordinate pairs P = {{xk,yk), {xk+iiVk+i)}) where P is 
a horizontal move if Xk+i = Xk and yu+i — yu + 1, a vertical move if x^+i = + 1 and 
£/fc+i = 2/fc, and a diagonal move if x^+i = + 1 and yu+i = yk + 1- Define a pat/i to as 
a sequence of coordinate pairs (x ,y ), (xi,yi), . . . such that (x ,yo) = (0,0) and each 
successive pair is either a horizontal move, a vertical move, or a diagonal move. The cost 
of a path is the product a h b d c v , where h, d, and v represent the total number of horizontal, 
diagonal, and vertical moves, respectively, in the path. Let f(i,j) represent the total cost of 
all paths to (i,j). Note that /(0,0) = 1, f{0,j) = a j for all j > 0, /(z,0) = c l for all i > 0, 
and moreover that if f(i, j — 1), f(i — — 1), and f(i — 1, j) are known, then the total cost 
f(i,j) can be computed as 

f(ij) = af(i,j - 1) + bf(i - 1, j - 1) + cf(i - 1, j). 

That is, M[i,j] = f(i,j) for all 

The total cost of all paths to can also be obtained as the sum, for k < i, of the cost of all 
paths to (i, j) that include exactly k diagonal moves, i — k vertical moves, and j — k horizontal 

3 



moves. For each k < i, there are 
j + i — k 



k 



ways to choose the columns for the diagonal moves, 



and then 



i — k 



ways to choose the locations of the vertical moves from among the 



remaining moves. Summing the cost a? k b k c l k over k < i yields the expression pi). □ 
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Throughout the remainder of the discussion below, we work with a fixed prime p, fixed 
integers a, b, and c, and a matrix M defined as in Definition 4A_ Since all arithmetic will be 
modulo p, we assume the entries of M are the residues modulo p. 



Definition 4.6 Let M(x, y, u) denote the finite uxu submatrix of M whose lower left corner 
is at (x, y); that is, 

M(x, y, u) [i, j] = M[x + i,y + j] 

for < i < u, < j < u. 

Definition 4.7 Let k > 0. A 1-block of size p k is the matrix M(0,0,p k ) . For < n < p, 
an n-block of size p k is a p k x p k matrix B for which B = nM(0, 0,p k ) mod p. 

Observation 4.8 Note that in the terms of the preceding definitions, M is numerically 
p- self- similar modulo p if and only if for all s,t < p and all k > 0, M \sp k , tp k , p k ) is a 
M[s,t]-block of size p k . Note also that an n-block satisfies modulo p, and therefore is 
completely determined by its first row and first column. 

Lemma 4.9 Let k > 0. Then M[0,p k - 1] = M[p k — 1,0] = 1 mod p. 

Proof. By definition, M[0,p k — 1] = a pk ~ l and M[p k — 1,0] = c pk ~ l . Using Fermat's little 
theorem and a simple induction on k, a v = a mod p, so aP _1 = 1 and similarly c p _1 = 1. 



Lemma 4.10 (a) Let k > 0, t < p, j < p k , andn = M[0,t] = a 1 . Then M(0,tp k ,p k )[0,j] = 



(b) Let k > 0, t < p, j < p k , and n = M[s, 0] = c s . Then M(sp k , 0,p k )[i, 0] = nd 
Proof. By definition 



where the latter is equivalent to a*a J using Fermat's little theorem as in the proof of the 
previous lemma. The second part is similar. □ 

Lemma 4.11 For all < i, j < p, 



□ 



na 3 



M(0,tp k ,p k )[0,j] = M[0 1 tp k +j] 




(a) aM[i - l,p- 1] +bM[i,p- 1] 



mod p. 



(b) bM\p-l,j-l]+cM\p-l,j} 



mod p. 
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Proof. Note that for < k < p, ( ? ) is divisible by p, so for i < p, 



k=0 



^p—kiki—k 



= a p c\ 



a p c l 



Then by definition, for < i < p, 

aM[i,p-l]+bM[i-l,p-l] + cM[i-l,p) = M[i,p), so 
oM[i,p - 1] + bM[i - l,p - 1] + c(a p c^ 1 ) = a p d and hence 
aM[i,p-l] + 6M[j-l,p-l] = 0. 

The argument for (b) is similar. □ 

Lemma 4.12 (a) Let < x < y and fix a row i > 0. Let n = M[i,x}. Suppose that for 
each column between x and y, the sum of adjacent entries of row i — when scaled by 
the coefficients b and c, is modulo p; that is, for all < j < y — x , 

bM[i - 1, x + j - 1] + cM[i -l,x + j}= 0. 

Then M[i, x + j] = na J for all < j < y — x 

(b ) Let < u < v and fix a column j > . Let n = M[u, j] . Suppose that for rows between 
u and v, the sum of adjacent entries of column j — I, when scaled by the coefficients a 
and b, is modulo p; that is, for all < i < v — u, 

aM[u + i,j - 1] + bM[u + i- l,j - 1] = 0. 

Then M[u + i, j] = nc % for all < i < v — u 

Proof, (a) For j = we have M[i, x+0] = n. Having shown for an induction that M[i, x+j] = 

M[i,x + j + l] = aM[i,x + j]+bM[i-l,x + j) + cM[i-l,x + j + l) (5) 
= aM[i,x + j] (6) 
= na J+1 . (7) 

The proof for (b) is similar. □ 



Proof of 4-2 We show the following by induction on k. Note that (a) implies that M is 



numerically p-self-similar. 
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(a) For s,t < p and n = M[s, t], 

M(sp k -\tp k -\p k - 1 ) = n- M(0,0,p k ~ 1 ) mod p. 
That is, M(sp k ~ 1 ,tp k ~ 1 ,p k ~ 1 ) is an n-block of size p k ~ x . 

(b) Let n < p and Let B be any n-block of size p k . Then for all < i, j < p , 

aB[i,p k - 1] +bB[i- l,p k - 1] = mod p, and 
bB\p k -l) + cB\p k = mod p. 

(c) Let n < p and Let S be any n-block of size p k . Then 

B[0,p k -1] =B[p k -1,0] = n modp 



£[p fc - l,p k - 1] = n mod p. (9) 



(d) Let s,t < p and let 



A = M((s-l)p k ,tp k ,p k ), 

B = M((s-l)p k ,(t-l)p k ,p k ), 

C = M((s- l)p k ,tp k ,p k ), and 

D = M(sp k ,tp k ,p k ). 

Suppose that A is an x-block, B is a y-block, and C is a ^-block and let w = ax+by+cz. 
Then D is a w-block. 

Note first that (d) follows from (b) and (c). That is, fix k and assume that (b) and (c) hold, 
and assume the hypothesis of (d). Then applying (c), 

A[0,p k -1] = x, (10) 
B\p k - l,p k - 1] = y, and (11) 
C\p k -1,0] = z, (12) 

and so £)[0,0] = w. By (b), the rightmost column of A satisfies the hypothesis of Lemma 
4.12 a), so D[i, 0] = wc % for i < p k . A similar argument shows that D[0, j] = wa^ for j < p k . 



Thus D[i,0) = wM[i,0) and D [0, j] = wM[0,j] for i,j < p k , so we have D[i,j] = wM[i,j] 
for all i, j, < p k by Observation |4.8[ hence D is a iy-block. 

We first establish the base step for k — 1. Part (a) asserts only that M[s, t] = M[s, t]-M[0, 0]. 



Part (b) follows from Lemma 4.11 For (c) ^is immediate from Lemma [4. 1Q| and [9] follows 



from (b) using the observation that for < i < p, 

M[i,p- 1] e 

(and if p — 2, p — 1 = 1). Part (d) follows in general from (b) and (c) as shown above. 



1 if i is even 

p — 1 if z is odd 
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Let k > 1 and assume the induction hypothesis holds. We first show that for t < p, 
M(0, tp h ,p k ) is an n-block of size p k , where n = M[0,t]. The t = case is the definition 
of a 1-block. Suppose we have shown that C = M(0, (t — l)p h ,p k ) is an m-block, where 



m = M[0,t - 1]. Let D = M(0,tp k ,p k ). We know from Lemma |4.10| that D[0,j] = no? 
for j < p k and in particular £)[0,0] = n. Since C is an m-block of size p k , by part (b) 
of the induction hypothesis the rightmost column of C satisfies the hypothesis of Lemma 



4.12[ b), and hence -D[i,0] = nc % for all i < p k . Since D[i, 0] = nM[i, 0] for i < p k and 



D 0,j] = nM[0,j] for j < p k , we have D[i,j] = nM[i,j] for all i,j < p k using Observation 
4.8 Thus D is an n-block. A similar argument shows that M(sp k ,0,p k ) is an M[s, 0]-block 



for each s <p k . 

The next step is to show, inducting on s and t, that M(sp k ,tp k ,p k ) is a M[s, t]-block. Let 

x = M[s,t-1], 
y = M[s-l,t-l], 
z = M[s - l,t), and 
w = M[s, t] = ax + by + cz. 

and suppose that 

A = M(sp k , (t - l)p k ,p k ) is an x-block, 

B = M((s-l)p k ,(t-l)p k ,p k ) is ay-block, and 

C = M{{s - l)p k ,tp k ,p k ) is a 2-block. 

Then by part (d) of the induction hypothesis, M \sp k , tp k , p k ) is a M[s, t]-block, which es- 
tablishes part (a). 

Next, we know from (a) that for each s < p, B = M(sp k , (p— l)p k ,p k ) is a M[s,p— l]-block. 
Using part (b) of the induction hypothesis, we have 

aB[i,p k - 1] + bB[i - l,p k - 1] = (13) 

for < % < p k . In addition, if s < p — 1 and A = M((s + l)p k , (p — l)p k ,p k ), then A is a 
M[s + l,p— l]-block, so by part (c) of the induction hypothesis, 

A[0, p k - 1] = M[s + 1, p - 1] and 
B\p k -l,p k -l] = M[s,p-1}. 

It follows that 

aA[0,p k - 1] + bB[p k - l,p k - 1] = 0. (14) 
Now consider any < i < p k+1 ; then 

aM[i,p k+1 - 1] + bM[i - l,p k+1 - 1] ee 0, (15) 
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which follows from (14) if i is a multiple of p k and from (13) otherwise. Then (b) is obtained 
by multiplying (15) by n. 



To establish (c) note that by Lemma |4T9 

M[0,p k+1 - 1] = M[p k+1 - 1,0] = 1. 

Also, B = M((p—l)p k , (p—l)p k ,p k ) is a M\p— l,p— l]-block by (a), where M[p— l,p— 1] = 1, 
so using part (c) of the induction hypothesis, 

M[p k+1 ,p k+1 ] = B[p k - l,p k - 1] = 1. 

Then (c) is obtained by multiplying the equivalences above by n. That (d) follows from (b) 
and (c) has already been shown above. □ 



We conclude this section with one example of a numerically self-similar fractal other than 
the Sierpinski carpet; in this case with a = 1, b = 2, c = 2, and p = 5. A portion of this 
structure is shown in Figure [4] with colors representing the 5 numerical values. 



5 Conclusion 



We have shown that the discrete Sierpinski carpet self-assembles in Winfree's Tile Assem- 
bly Model and, moreover, that it is an instance of an infinite class of discrete fractals that 



self-assemble. The key ingredient of this result was Theorem 4.2 that certain recursively gen- 



erated infinite matrices have a strong self-similarity property, which we defined as numerical 



self-similarity. Theorem 4.2 is a strong generalization of known results on self-similarity in 
Pascal's triangle which underlie the previous work on self-assembly of the discrete Sierpinski 
triangle. 

The recursively generated matrices we study provide a rich source of examples of self-similar 
fractals, despite the obvious simplicity of the linear function used to generate them. We 
have investigated matrices generated using more complex recursive relationships and the 
preliminary results are inconclusive; that is, although there are examples that appear to 
generate fractal structures (subsets of the plane with dimension strictly less than 2), none 
of the structures observed so far is self-similar. 

The discrete fractals we have investigated all self-assemble in a "progressive" way; that is, 
a tile that binds at location (x, y) is always determined by tiles at locations (x',y') with 
x' < x and y' < y. An interesting question is that of finding an exact characterization of the 
self-similar fractals that can be tiled progressively. A restriction to progressive tiling rules 
out self-similar fractals with blocks of zeros along either axis. It remains open whether every 
symmetric, self-similar discrete fractal without zeros along the axes can be tiled progressively. 
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